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1 Introduction 

Magnetars are fascinating objects that are thought to be neutron stars powered by 
their strong internal magnetic fields [Tj. 

There is evidence that neutron stars suffer a long term spin-down. Moreover, 
many sudden spin-ups, known as glitches in the literature, have been observed in 
pulsars and magnetars M. 

Clear evidence of a sudden spin-down was detected in the Anomalous X-ray Pulsar 
AXP IE 2259+586, an object cataloged as a magnetar [TJ. This event received the 
name “anti-glitch”. To adjust the timing data from Swift, two different interpretations 
for the observational evidence were proposed: (i) an anti-glitch in which Xv/v = 
—3.1(4) x 10 - ' followed by a spin-up event of amplitude Xv/v = 2.6(5) x 10 7 ; (ii) 
an anti-glitch in which Xv/v = —6.3(7) x 10 -7 followed by a second anti-glitch in 
which Xv/v = —4.8(5) x ICC' [4j. Based on a bayesian analysis, model (ii) is favored 
0 . 

Regarding the energetics related with this event we mention that consistently with 
the epoch of the anti-glitch, Fermi/G BM detected a hard X-ray burst with a duration 
of 36 ms [6]. The observed fluence in the 10-1000 keV band corresponds to an energy 
release of E 7 ~ 10 38 erg. Moreover, an increase by a factor 2 in the 2-10 keV flux 
was also observed [4], resulting in a Ex ~ 10 41 erg energy release [7]. 

AXP IE 2259+586 has a characteristic age of 10 6 yr, a ~7 s period and a spin- 
inferred surface dipolar magnetic held of B^ ~ 5.9 x 10 13 G0. Over the last two 
decades it has been monitored by the Rossi X-ray Timing Explorer and the Swift 
X-ray Telescope. With the exception of two spin-up glitches in 2002 [8] and 2007 [9], 

1 http://www. physics, mcgill.ca/^pulsar/magnetar/main.html 
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a timing event in 2009 [9] and this anti-glitch in 2012, the source showed a stable 
spin-down rate. 

Several explanations for this anti-glitch event have been proposed, which can be 
divided in two different families of models: one based on an external |Zl lEH HD 02] 
origin and the other on an internal HU H3] one. Despite several searches in radio 
and X-ray wavelengths, no surrounding afterglow was detected [4], arguing against 
a sudden particle outflow or wind-driven scenario. In this sense, a more promising 
approach in order to explain the phenomenum seems to be an internal rearrangement 
of the star. 

In this work we present a simple internal mechanism which could account for the 
observed sudden spin-down of the star |14j . The central idea behind this model is 
that as a consequence of the natural long term decay of the internal magnetic field, 
an initially prolate-shaped stable stellar configuration becomes unstable enough to 
crack the crystallized stellar crust. Then, the re-accommodation of the star into a 
stable more-spherical shape, could naturally lead to the occurrence of an anti-glitch, 
as a consequence of the conservation of angular momentum. A similar scenario was 
also suggested to account for the SGR 1900+14 event [ T5] . 


2 The proposed mechanism 


Deformations for a rotating uniform-density with a mixed poloiclal-toroidal magnetic 
field configuration were calculated in [T6l l?j. Studies with a wider and more realistic 
family of equations of state were performed recently by |T8] . As in both set of works, 
the quadrupolar distortions obtained are of the same order of magnitude, in our 
model, we used a simple uniform density star which allows to perform analytical 
calculations. In addition, since magnetars like AXP IE 2259+586 show long spin 
periods, deformations due to rotation are negligible, and thus, we do not consider 
them on our model. All these different works conclude that equilibrium configurations 
for stars with strong poloidal magnetic fields or in rapid rotation, are oblate while, 
when internal toroidal fields dominate the magnetic field configuration, prolate stars 
result favored. 

For incompressible stars of uniform density and mixed poloiclal-toroidal magnetic 
field, the quadrupolar distortion of equilibrium configurations, e, of the l = 2 volume 
preserving mode is given by: 


e = 



25 R 4 (, 2 v 
24 G n M 2 V t] 



( 1 ) 


where M and R are the mass and the radius of the undeformed star, respectively. Gn 
is the gravitational constant and (B^ p ) is the mean value of the square of the toroidal 
and poloidal magnetic field strengths. 
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A purely toroidal magnetic field configuration is known to be unstable [19], but a 
poloidal component with energy E p /E t = B 2 /B 2 ~ 1 — 5% can stabilize the magnetic 
held configuration [20]. Thus, we neglect the poloidal contribution to equation (Q]) 
considering that ( B 2 ) <C ( B 2 ). 

Neutron star crusts are strong enough to support cllipticities up to a critical value 
of e c < 4 x 10~ 6 before cracking [2Tj . 

After these general considerations we present the theoretical picture that want to 
explore: given an “initial”, mostly toroidal, magnetic held with strength ( B \), the 
neutron star crust crystallizes in a prolate equilibrium configuration, with an elliptic- 
ity e~ given by equation (JI]). A series of effects that are thought to take place inside 
the neutron star produce a progressive decay of the original magnetic held in time 
scales ~ 10 5 yr [22]|23j. As a consequence, the prolate configuration with e~ departs 
from equilibrium until the stellar crust reaches a critical strain and cracks. Then, the 
stellar structure achieves a new stable and less prolate configuration, with cllipticity 
e + , associated to the present or “hnal” magnetic held strength, (B{). Because the 
more spherical configuration, e + , has greater moment of inertia respect to the spin 
axis, chosen to be z in our case, with respect to the previous e~ configuration, and 
considering that in the absence of an external torque angular momentum conserves, 
this sudden change in the stellar structure can easily account for the observed sudden 
frequency spin-down. 

A change in the oblateness of a uniform density star induces a spin frequency shift 
given by 


Av (1-2 e+1 1 / 3 2 

-= ----1 « - e~ — e 

v (l-2e-)V3 3 6 


( 2 ) 


relationship that, as a function of (B^) 2 , can be written as (for details on the calcu¬ 
lations see m 


Av 

v 

from where it follows that 


2 25i? 4 

3 24 G n M 2 




(3) 


(B{) » \J(Bi) 2 - Bl (4) 

where B q = — | ^ 24 ^ 4 f would certainly change for a deeper treatment of the stellar 
structure. 

In this sense, it would be interesting to perform more detailed calculations con¬ 
sidering, for instance, a neutron star composed by a solid crust surrounding a liquid 
core, as the one performed to study starquakes in rotation-powered pulsars [24]. How¬ 
ever, at this point, we prefer to keep on the simplest model to present here order of 
magnitude calculations to give confidence to the suggested scenario. 
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Figure 1: Physical solutions to equation (J4]) for Ai//z/ = —6.3 x 10“' as function of 
( B \), for three different neutron star configurations (see the legend). In black we plot 
the identity function as a reference. 


To estimate the energy released by this mechanism during the anti-glitch event, we 
use the classical model developed for pulsar glitches [25]. Three energy contributions 
are considered: gravitational, coming from the global change in the stellar shape; 
rotational, associated to the spin frequency shift, Az//z/, and strain tensions released 
by the crust. This strain energy is accumulated because, despite that the magnetic 
held decays from ( B l t ) 2 to (B{) 2 , the crystallized crust keeps its original shape of 
by increasing its internal tension, departing from equilibrium. Once the critical 
strain is achieved, the crust cracks and the star re-accommodates into an e + equilib¬ 
rium configuration, releasing the stored energy. As the length scale, £, associated to 
the change in shape in our model is small l/R ~ e, we assume that magnetic held 
reconnection does not take place. For this reason, even though we are treating with 
a magnetar, we do not consider the magnetic energy release due to the displacement 
of magnetic held footpoints, following the standard Soft Gamma Repeater picture 

mm- 


3 Results 

In the framework of our model, we estimate the long-term decay in the magnetic 
held strength needed to account for the [S.vjv observed in the anti-glitch of magnetar 
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|A<B t >|/<B'> 


Figure 2: Avjv values obtained from equation ([3]) as a function of the relative 
change \A(B t )\ /(B l t ) in the mean toroidal magnetic held strength. As a refer¬ 
ence we also show the Av/v values of the anti-glitch/anti-glitch pair observed in 
AXP IE 2259+586 in shaded rectangles with their corresponing error bars. 

AXP IE 2259+586. 

As a function of the mean initial toroidal magnetic held strength, (B l t ), and 
for three different neutron star configurations, we plot in Figure ([I]), the physi¬ 
cal solutions, A(B t ) < 0, to equation (]4]). For this, we assume a frequency jump 
Avjv = —6.3 x ICE' equal to the hrst of the two events of model (ii). In black we 
plot the identity as a reference. 

For a typical neutron star with a mean toroidal magnetic held of ( B t ) = 2 x 

10 15 G, which corresponds to a maximum value for the magnetic held strength Bm > 

10 16 G [20], we estimate that a decay in the magnetic held of about ~10% could be 
responsible for the observed spin-down. These qualitative result is almost insensitive 
to other acceptable values for mass, radius and magnetic held strength of magnetars 
(see Figured]). Detailed studies of the magnetic held evolution in neutron stars show 
that a magnetic held decay of ~10% is easily achieved after t < 10 6 yr for a magnetar 
like AXP IE 2259+586 [23]. 

For the adopted neutron star configurations, we find that a minimum value for 
( B l t ) (different in each case) is needed in order to have a solution to equation ([3]) for the 
observed Avjv. This critical value is, in any case, several times 10 14 G, which avoids 
the occurrence of anti-glitches in normal pulsars, only allowing this phenomena to 
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occur in strongly magnetized neutron stars, i.e. magnetars. This result can be used 
to explain why, despite many pulsars have been thoroughly monitored for several 
decades, no sudden spin-down event of this kind has been detected at all. 

In Figure (EJ we present Av/v from equation (J4j) as a function of the change in the 
mean toroidal magnetic held strength. The shaded rectangles represent the values of 
A v/v from model (ii) and their corresponding error bars. 

We showed how our simple model can be used to explain the timing behaviour of 
the anti-glitch observed in magnetar AXP IE 2259+586. Hence, we now focus our 
attention to analyze the energetics. 

For a “typical” neutron star and assuming a normal neutron star crust, we estimate 
a gravitational energy release of ~ 10 42 erg, while the rotational and crustal strain 
contributions are several orders of magnitude smaller. Assuming a typical SGR, the 
magnetic energy release is ~ 10 41 erg. Hence the total energy emitted after the anti¬ 
glitch in gravitational waves, particles and electromagnetic radiation should be of the 
order of 10 42 erg which is compatible with the observations. 


4 Conclusions 

We present here a very simple model based on an internal mechanism to explain the 
anti-glitch observed in magnetar AXP IE 2259+586. We suggest that a long term 
magnetic field decay of approximately ~10%, from an initial ( B\) ~ 10 15 G, would be 
enough to de-stabilize an originally prolate stellar configuration cracking the neutron 
star crust to a “more spherical” one. As a result, the sudden change in the moments of 
inertia of the star would naturally lead to a sudden spin-down, as the one observed in 
|4|. Under this scenario, in addition, considering a “typical” neutron star, we estimate 
an energy release of ~10 42 erg, which is in agreement with the emission detected by 
Fermi and Swift observatories in the epoch when the anti-glitch occurred. 

As a corolary, our simple model predicts that an anti-glitch as the one detected 
by m in AXP IE 2259+586, can only be achieved if the mean internal toroidal mag¬ 
netic field of the neutron star is several times 10 14 G, as in magnetars, avoiding the 
occurrence of anti-glitches of this amplitude in normal pulsars, which might be the 
reason why this event is the first of its kind to be detected. 
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